NATIONAL ADVISORY COMMITTEE 
FOR AERONAUTICS 

TECHNICAL NOTE 3594 


EFFECT OF TRANSVERSE BODY FORCE ON CHANNEL FLOW 
WITH SMALL HEAT ADDITION 
By Simon Ostrach and Franklin K. Moore 


Lewis Flight Propulsion Laboratory 
Cleveland, Ohio 




3593 


TECH LIBRARY KAFB, NM 



RATIONAL ADVISORY COMMITTEE FOR AERONAUTICS 
TECHNICAL NOTE 3594 

EFFECT OF TRANSVERSE BODY FORCE ON CHANNEL 
FLOW WITH SMALL HEAT ADDITION 
By Simon Ostrach and Franklin K. Moore 


SUMMARY 

The steady, compressible, inviscid channel flow to which heat is 
added at a cross-sectional plane and which is also subject to a trans- 
verse body force is analyzed. The prime parameters governing the flow 
. are a dimensionless heat parameter and the Mach and Froude numbers. So- 

H lutions which are qualitatively similar are obtained for the cases in 

N which either the Mach number is small or the Froude number is large. If 

the heat is added uniformly over the plane, the streamlines are dis- 
placed in the direction of the body force, and downstream of the heat 
addition the flow is found to be a shear flow. If the heat addition is 
concentrated nqar the center of the cross-sectional plane, a strong vor- 
tex motion appears downstream in addition to the shear flow. In the 
first part of the analysis, the through- flow Mach number is considered 
small. Subsequently, the alternate assumption of large Froude number 
is made . For large Froude numbers it is found that an increase in Mach 
number to unity increases the magnitude of the body- force effects. 

The configuration studied can be considered to be an idealization 
of, for example, the flow in a ram jet mounted at the tip of a whirling 
helicopter rotor blade. The effects of the transverse body force are 
important in this example. Actually, even if the Froude number is 
large, the vortex motion due to any nonuniformity of heat addition may 
be quite vigorous . 

INTRODUCTION 

It is well known that the motion of a fluid in which there are den- 
sity variations is altered by the action of body forces on the fluid. 
Until recently, such buoyancy effects have been negligible in practical 
applications, because both the density variations and the body forces 
we re relatively small. However , considerations of present-day propul- 
sion systems with their associated large heat additions, rapidly rotat- 
ing components, or their operation in sudden maneuvers at high speeds 
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suggest a reevaluation of these effects. Of prime importance in this 
respect is the effect of the body forces on the internal aerodynamics 
of such propulsion systems or their components . The influence of the 
body forces has been studied for channel flows with uniform heat addi- 
tion and with the body force para ] 1 el to the direction of motion (see 
ref. 1, e.g.). The flow with heat addition in pipe bends , through 
rotor-tip combustors of helicopters , or through ram jets in abrupt ma- 
neuvers represents configurations of interest in which, on the other 
hand, the body force is transverse to the flow direction. 

Configurations of the latter type will be considered herein. More 
specifically, in the present analysis a steady uniform (in velocity and 
temperature) flow in a cha nn el subject to a transverse body force is 
assumed perturbed by a steady small addition of heat at a cross- 
sectional plane (fig- l) . The fluid is considered to be inviscid. The 
neglect of viscosity is in contrast to the analysis of reference 1. The 
reason for this difference can be explained as follows: 

Important body- force effects appear in regions of large temperature 
gradients transverse to the body- force direction. In internal natural- 
convection flow problems (ref. 1, e.g.), the largest temperature gradi- 
ents occur in the viscous heat-conducting layers of fluid next to the 
channel walls; accordingly, a parameter (Grashof number) comparing vis- 
cous and body- force effects governs the solutions of these problems . 

In the present problem, the important temperature gradient is not asso- 
ciated with viscous wall effects, but, rather, with the addition of heat 
at a flame front or combustion zone. Therefore, the fluid may be con- 
sidered inviscid in this problem; and the Froude number, which compares 
inertia forces and body forces, and not the Grashof - number is of chief 
importance. In addition, of course, a parameter specifying rate of heat 
addition to the fluid and the Mach number are also of significance. 

If heat is added uniformly over the cross section (fig. 2(a)), the 
sudden reduction of density at x = 0 (for small through Mach number) 
would be expected to cause the streamlines to have a tendency to rise 
in the vicinity of the heat addition. Perhaps the most interesting ef- 
fects may be expected when heat is added over only part of the cross 
section, as in figure 2(b) . Then, the difference in buoyancy between 
the flow near the center and the flow near the edges of the channel 
would be expected to result in the formation of vortices with axes 
alined with the through flow. 

With a view to studying these flow patterns, solutions are sought 
for the perturbation equations further simplified by the assumption of 
sTna l 1 Mach number or, alternatively, of large Froude number. 
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DIFFERENTIAL EQUATIONS AND BOUNDARY CONDITIONS 
Basic Equations 

s'" 

The differential equations governing the steady flow witii heat ad- 
dition of a compressible, inviscid, nonconducting fluid which is subject 
to a body force are, in rectangular Cartesian tensor notation, 


o 

II 

? 

(la) 

TT Su i _ SP 

pu j 3xJ “ pf i ' 5x7 

(lb) 

St du-j 

pc v u j Q " P 3xJ 

(lc) 

- = RT 
P 

(Id) 


(A complete list of the symbols used herein is presented in the appen- 
dix.) Equations (la) to (lc) express, respectively, the conservation 
of mass, momentum, and energy; equation (id) is the state equation. 

The more familiar notation U,V,W and X,Y,Z instead of 

and X-pXgjX-j will be used in the remainder of this report. 


Boundary Conditions 

Particular consideration is given to the configuration illustrated 
in figure 1: Far upstream, there is a flow with uniform velocity and 

temperature in a straight channel, so that 


U(-®,Y,Z) =* U w (a constant) 

(2a) 

V(-°°,Y,Z) = W(— , Y,Z) = 0 

(2b) 

T(-®,Y,z) a 0^ (a constant) 

' (2c) 
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The channel has a uniform rectangular cross section, so that 

V(X,0,Z) = V(X,H,Z) = 0 (3a) 

W(X,Y,0) =» W(X,Y,L) = 0 (3b) 

The fluid is subject to a transverse body force f2 acting in the 
direction of positive Y. Thus, f]_ and are defined to be zero, 

and ±2 is a negative quantity if the body force acts in the direction 

illustrated in figure 1. The action of the body force leads to pressure 
and density gradients in an equilibrium state. This equilibrium will be 
disturbed by the addition of heat to the fluid. 

Thus two specific problems are treated herein according to the 
specified distribution of heat addition. If heat is added uniformly 
(fig. 2(a)) over the cross section at X = 0, the problem is two di- 
mensional in X and Y. A three-dimensional problem will also be 
treated, in which the rate of heat addition at X = 0 varies with 
the coordinate Z, which is transverse to the direction of body force 
(see fig. 2(b)). 

However, the analysis will be carried forward as far as possible 
with an arbitrary function to represent heat addition. 


Linearization of Equations 


In the present analysis, a perturbation procedure will be employed. 
Accordingly, the actual rate of heat addition per unit volume Q at 
X = 0 is included in the dimensionless quantity 


eq(X,Y,Z) 


QH. 


p U c„T 

oo cd p a 


( 4 ) 


where the function q is of unit order and e is a small parameter. 

As a consequence of the ama.1 1 rate of heat addition, the flow is assumed 
to be disturbed only slightly from the uniform conditions prevailing far 
upstream (eqs. (2)). The disturbances are assumed to be of the order 
e, and the following definitions are made: 



T = T ro (l + ex), P = + eP r n, p = + ep r X 


( 5 ) 
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where p r is a constant reference density. Also, dim ensionless coordi- 
nates sire defined as 

x,y,z = X/H, Y/H, Z/H (6) 


respectively. 


Substituting equations (5) and (6) into equations (l) yields, to 
zero order in e. 


so that 


dp yM 2 

oo ao _ 

dy“ ~ " ^2 ® 

p, (y) = KEP. (y) 



( 7 ) 


p = p exp(-yM 2 y/Pr 2 ) 

CO X OO 


and 

P ot = P r exp(-yM^y/Pr 2 ) 


These expressions involve the Froude and Mach numbers defined by 

Pr = U > /V(-f 2 )=” 

► 

M. 5 


( 8 ) 


Clearly, equations (7) egress the conditions of static equilibrium of 
density and pressure with body force, if no heat is added. In such a 
case, the presence of a uniform through velocity would not affect 

this equilibrium . 


The following equations to the first order in e result upon elim- 
ination of u, t, and X from the equations obtained by combining (l), 
(5 ) and (6 ) : 


^XZ “ V XX = 0 

(9a) 

u xy “ v xx + X ( u x + v y + v z) = 0 

(9b) 

(1 - M 2 )ux + v y + w z + \(y - l)M^v = q 

(9c) 
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■where X = l/Fr 2 and, hence, represents essentially the ratio of "body 
to inertia forces. 


Definition of a Potential 


In order to reduce the number of variables of equations (9), a 
function ijr is defined such that 

*xx = + v y + v z) ( i0 ) 

From equations (la), (5), and (10) it may be shown that the func- 
tion ljr is related to the perturbation in density: 


% 


-XX 


Then, in view of equations (2) and (10), the Cauchy-Riemann conditions 
may be extracted from equations (9) in the form 


^ " w x 


Uy - (v - *) x = 0 ^ 

Wy - (v - I{r) z = 0 

Accordingly, a potential cp may be defined so that 

u = q, x 

v - ijr = (p 5 

v= cp z 

Equations (9c) and (lO) to (12) yield 


(H) 


( 12 ) 


V 2 ^ - M 2 ^ + (r - l)XM 2 cp v = q. - tr v - (r - l)XM^r (l3a) 

(13b) 


^xx ' ^y = Xv2( P 


Equations (2), (5), (6), (10), and (12) yield the boundary condi- 
tions to which equations (l3) are subject: 

cp(-co,y,z) = t(-®>y>z) = 0 


(14a) 
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<p y (x,0,z) + i|r(x,0,z) = <p y (x,l,z) + t(x,l;z) = 0 (l4b) 

<P z (x,y,0) = <P z ^y,|^= 0 (14c) 

Hereinafter, the heat addition will he assumed to occur abruptly at 
x = 0. That is. 


q S 6(x)q*(y,z) 

where 6(x) is the Dirac delta function. 


(15) 


SOLUTIONS OF SPECIAL PROBLEMS 
Small Mach Number 

Equations (13), while linear, are not especially simple. There- 
fore, results will be obtained under further simplifying assumptions 
If the Mach number is sufficiently sma1 1 so that « 1 and the 

CO 

Froude number is of unit order, then, in place of equations (13), 

P x o x 

i|r » x I d x / 5(x)q*dx - Ax l(x)q* (16a) 

\J-co co 

= 6(x)q* - Ax l(x)q^ (l6b) 

where l(x) is the unit step function. 

Quite evidently, equations (13) may also be simplified by the as- 
sumption of small (but nonvanishing) A.. This approach will be discussed 
in the subsequent section. Large Froude Number. 

Reduction of potential problem to familiar forms . - In order to ob- 
tain a more familiar problem of potential theory, equation (16b) may be 
solved in two parts, defining 


<P = <Pq + Tg 


(17a) 


V 2 ^ = q = 5(x)q* 


where 


(17b) 
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and 



-Xx l(x)q* 


(17c) 


A function which, is a particular integral of equation (17c) for 
x > 0 is 


<P 3 = -Xx l(x)h(y,z) (l8a) 

where h(y,z) satisfies Poisson's equation 

kyy + h zz = <ly ( 18b ) 

hut is otherwise unspecified, as yet. Then, if the definition 

^2 = ^3 + ^4 C 19 ) 

is made, substituting equations (18) into equation (17c) yields 

V 2 <P 4 = XS(x)h(y,z) (20) 


Except precisely at x = 0, (p^ and q> 4 satisfy Laplace's equa- 
tion. Solutions of Laplace's equation, valid for x < 0 and x > 0, 
are denoted by superscripts - and +, respectively. Integrating 
equations (17b) and (20) across the discontinuity in the plane x => 0 
yields the compatibility conditions 


("k ■ "yLo 

(p4x " (p 4j) xss 0 


= q* 

(21a) 

= Xh 

(21b) 


In addition to equations (2l), the boundary conditions (eqs. (14)) must 
be satisfied. 


Thus, after the singular solution (eq. (18a)) is extracted, the re- 
maining differential equations are . simply Laplace ' s equation (applied 
independently for x > 0 and x < 0) subject to the jump conditions of 
equations (2l) and Poisson's equation (18b) in two dimensions. 
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Fourier series solution for g* = q~*(z) . - The foregoing equations 
may readily be solved by Fourier analysis, if q* is assumed to be a 
function only of z (fig. 2), representable as 


q*(z) 


2 

n?=0 


^cos jmaz 


(22) 


where a « H/L . A solution of Laplace ' s equation which also satisfies 
equations (l4a), (14c), and (21a) is 



V 


CO 

-z 

n=l 


On 

2jtna 


e -jtnax cog jmaz 


(23a) 


w 

i 

t>3 

O 


<P 


1 



n=l 


gimaxQQg 


(23b) 


Because 


q* => 0, it is permissible to take 



sinh jtna 



jtna cosh 


jtna 


cos Jtnaz 


(24) 


as the proper solution of equation (18b) . In view of the definition of 
ijr (eq. (16a)) and of equation (22), this choice provides that q> 3 (eq. 

(18a)) satisfies equation (14b). Also, equation (14a) remains satisfied. 


The function cp^ must be found so that it is consistent with equa- 
tion (21b). Also, of course, cp^. must satisfy equations (14a) and (l4c) 
and together with (p^ must satisfy equation (14b) . 

Being antisymmetric about the plane y = l/2, the y-dependence of 
h(y,z) in equation (24) is represented by the Fourier series in order 
that the boundary conditions can be satisfied. Thus, 


sinh jtna 


Qr 


M) 


jma. cosh 


jtna 


•co 

z 

m=0 


B^cos mjty 


(25a) 
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where , by Fourier's theorem, 

= 0 (m zero or even) 
4 


~2T2 — rn ^ odd) 


1 


Then, equation (24) becomes 


(25b) 


CO CO 


h = y? y? BrrmCOS mrty cos itnaz 
ntaO n=0 


and, from equation (21b), the proper equation for q>. must be 


(26) 


it 

*4 


CO oo 



m=0 n=0 


®mn' 


Tux. 


4 


m 2_2 2 
m +a n 



a 2 n 2 


cos nuty cos imaz 


(27) 


satisfying equation (21b). 

To recapitulate equations (17a), (l8a), (19), (23), (24), and (27), 
the result for cp is 


X l(x) 


at 


rt=«0 

+ a£a 2 


sixth ttna 


tJl 


cos irnaz 


itna. cosh 


2«na 


^ e : ' : ’ tnct3<: coB jtnaz 


n=l 


jj EL - ■ e Trtx 'V^+“^cos Duty cos itnaz 


m=0 rfc*0 


( 28 ) 


Result for uniform heat addition . - In the instance of uniform heat 
addition (fig. 2(a)), q*(z) = Qq from equation (22), with all other Q n 

being zero. Then, incorporating equation (25b) and using L'Hospital’s 
rule yield for equation (28): 



l(x) 




e ' : F(2j+l)rtX 

2^ (2j + l) 3 

3=0 


cos(2j + l)ity> 


(29a) 
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Also, from equation (16a), 


i|r =» Xx 1 (x)Qq 


(29b) 


The velocity components, from equations (12), are 


O) 

tffca * V 

|i(x) 

1 - 4 - 1) 

Ol 

CD 

W 


L 




e T(2o+l)jtx 
(2j + l) 2 


cos(2j + l)jty> 


(30a) 


v^ 



e T(2j+l)nx 

— -o sin(2j + l)ity 

(2J + l) 2 


(30b) 


The result for vertical velocity (eq. (30b)) indicates that fluid 
streamlines are displaced downward as they traverse the section where 
heat is added in proportion to the amount of heat added Qq and the 

Froude number (through X) . Such a streamline is shown in side view 
in figure 3(a) . Inasmuch as equation (30b) is symmetrical in x, one- 
half the displacement occurs upstream of the heat addition. The total 
displacement of a streamline near the center of the channel (y = l/2) 
is approximately, from equations ( 4 ), (5), (7), and (30b), 



1 

8 


4-co 

X 


QH dx 
P U cmT 

00 00 jr 1 00 


(31) 


where the integral is, essentially, the total rate of heat addition per 
unit area. 


This displacement occurs in a distance of the order of H and is 
in a downward sense.- Figure 3(a) shows that downward displacement is, 
in fact, consistent wL th a buoyant force acting at x = 0. The stream- 
line curvature which is negative for x < 0 and positive for x > 0 
suffers an abrupt upward increase at x = 0- Of course, the restraint 
imposed by the walls at y =» 0,1 determ in es the particular streamline 
pattern to accompany the upward buoyant acceleration at x = 0. 

The result for the dimensionless through velocity u/tJ , from equa- 
tions (5) and (30a), is indicated in figure 3(b). Far upstream (station 
l) the flow has a dimensionless uniform velocity 1. Far downstream 
(station 4), the average dimensionless velocity is higher by eQ q 
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because of the heat addition.- 1 - Because of the buoyancy effect, the ve- 
locity profile has a uniform shear of such a magnitude that the dimen- 
sionless velocity difference between the bottom and top of the channel 
is cQqX- For large distances downstream, this is the dominant effect 

of Froude number. The sense of the shear is such that the through flow 
is more rapid in the constricted flow in the lower part of the channel 
(cf. figs. 3(a) and (b)), and less in the expanded flow at the top. 

This result is in accord with the usual behavior of one- di mensional low- 
speed flow in a channel of varying area. 

In view of equation (21a), the actual velocity jump at x = 0 must 
be independent of y and z. This requirement is met by the potential 
flow described by the series terms of equations (30) . Evaluated at 
x ss 0, these terms yield the dashed-line profiles in figure 3, which are 
parallel and differ by the constant bQq, between stations 2 and 3. 

The potential-flow terms attenuate exponentially in distances of 
the order of H, both upstream and downstream from the zone of heat ad- 
dition. At large distances, the most persistent of the series terms is 
the first one, yielding contributions to u and v proportional to 
exp [Tjtx] [cos ity] and exp [^x] [sin «y], respectively. 

Result for cosine-wave heat addition . - For the problem of cosine- 
wave heat addition (fig- 2(b)), 

q*(z) = — Q(l - cos 2itaz) (32) 

Cm 

Thus, from equation (22), Qq = ^ Q, Qg = - Q, and, from equation (28), 


|X l(x) 

1 1 

sinh 2ita^y - ^ 

cos 2?taz^ 

+ ■ 2jtaz 

1 " M y " 2 

2jta cosh jta 

4 ita 

L 

\ 


— 

— 


'sT' 5m0 grfttmx^g _ A cos 2«az — = 

m=tO - 11130 v 


®m2 


e : Boc'^4a 2 -Hn 2 cos 


+ (2 a) £ 


(33a) 


• ^The factor p /p is omitted here and in similar expressions in 
r 7 co 

this section on sma,! 1 Mach numbers, because its value is unity to con- 
sistent order in Mach number (see eqs. immediately following eq. (7)). 
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and, from equation (16a), 

IS/ 

i}r = Xx l(x) — (l - cos 2jtoz) (33b) 

Again, part of the solution involves potential functions, attenu- 
ating at distances of the order of H and L upstream and downstream 
of the station x = 0, where heat is added. 


On an average across the channel, streamlines are displaced down- 
ward, as in the previous case of constant heat addition. The velocity 
components -which dominate for larger x are, from equations (33) and 
( 12 ), 


u = \ Cftl - X 


sinh 2jta 


y - 2 " 


( y - 1) 


v = - g Q Xx 


2ita cosh ita. 
cosh 2rta/ 


cos 2itaz 


1 - 


( y - t) 


w = - i Q Xx 


cosh ita, 
sinh 2jta( 


cos 2itaz 


( y - s) 


cosh ita 


sin 2rtctz 


(34a) 


(34b) 


(34c) 


Equations (34) represent an incompressible vortex motion, satisfying the 
simple continuity equation u x + v y + w z = 0, and having the dimension- 
less components of vortlcity 


<2 X = v y 


v„ = -jtQxXa sin 2itaz 


2 y = u z - v x 


S z =v x -Uy = ^(l- cos 2naz) 


0 


(35a) 

(35b) 

(35c) 


Equation (34a) defines the part of the profile of u (fig. 4) which 
corresponds to the shear flow (eq. (3l)) occurring in the case of uni- 
form heat addition (station 4, fig. 3). This portion of the profile 
does not vary with x. In fact, from equations (35b) and (35c) the re- 
lated vorticity components S2 y and 2 Z are also constant with x. Of 

course, equations (34b), (34c),- and (35a) are valid only for moderate 
x, inasmuch as the present analysis is restricted to small perturbations. 
Since, physically, the vortices, must be finite even far downstream, the 
conditions far from the plane of heat addition must be determined from a 
more complete analysis. 
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In the present three-dimensional case, there is a tendency to form 
strong vortices with axes alined with the flow. These vortices acre in- 
dicated hy the streamline pattern for a cross section shown in figure 5 . 
The axial vorticity 2^. increases linearly with x, as equation (35a) 

evidences; the amount of circulation in each half of the cross section 
is found to he equal to eXQU^Hx. 

Another view of the vortex motion is obtained by tracing the vortex 
filaments of the flow. Each filament lies in a plane of constant y 
(2y = 0, from eq. (35b)) and is horseshoe-shaped, with the open end 

downstream. The vortex filaments lying in the plane y = l/2 are shown 
as dotted lines in figure 5. If the horseshoes are spaced equally in x 
(this is required by the fact that <2 Z is independent of x, from eq. 

(35c)), the number of filaments passing through a cross section appar- 
ently is directly proportional to x. Accor ding ly, 2 X must increase 

O 

linearly with x. These results are, of course, for essentially 

equal to zero. The effects of Mach number will be determined in the 
subsequent section. 


Large Froude Number 


As has been stated, simplification of equations (13) can also be 

obtained by assuming that the Froude number is large so that X = ~~ 

Fr d 

rather than the Mach number is srnal 1 relative to unity. Therefore, 
expanding q> and i]r in a Maclaurin series in X as 


q> = cp^ + X<p ^ + 

l{r at Xt^ + 

yields to zero order in X 

V 2 <p(°) - = q 

co XX. 



(36) 


(37) 


^The vortex system Illu strated in fig. 5 suggests a comparison with 
the system of bound vortices of a low-aspect-ratio wing, with a leading 
edge at x = 0 and with a span of L. The sense of the vortices in the 
present problem corresponds to a wing with negative lift distribution, 
and this is consistent with the downward inc lin ation of the flow which 
has previously been discussed as occurring in the vicinity of the heat 
addition. 


3593 



3593 


MCA TN 3594 


15 


o 


— 1 |/°^ = 0 
XX 


( 38 ) 


and to first order in X 

vV 1 ) - I^,W = - jyt 1 ) + ( r . l) M 2^p(°) + ,(0| 

+« . V 2 ,,^ + *<°> 

XX y 

The boundary conditions for equations (37) to (40) are 

<P^(-“,y,z) = (-oo^y^z) = 0 

Ty^CxjOjz) + i|r^(x,0,z) sip^(x,l,z) + | (l) (x,l,z) 

^^(x^O) = t f > z 1 ^( x ^ y ^l) = 0 


(39) 

(40) 

(41a) 
0 (41b) 

(41c) 


where i is 0 and 1 for the zero- and first-order equations, respec- 
tively. In view of equations (36), the zero-order solutions describe 
the flow with heat addition and no body force, and the first-order so- 
lutions represent the contributions due to the action of the body 
forces. Thus, in order to find the effects of a transverse body force 
on the flow, the first-order solutions must be determined. However, 
in the case of the small Mach number previously treated, the body- force 
effects are evident from the zero-order (in Mach number) solutions. 

Thus , for small Mach numbers , the body-force effects can be of the same 
order of magnitude as the effects of heat addition; whereas, for rhib.! 1 
X, the body-force effects must be small relative to the heat-addition 
effects. The solutions of the previous section are, therefore, perhaps 
more meaningful, but the present section serves to show how the Mach 
number affects the flow. 


Fourier series solution . - From the boundary conditions and equa- 
tion (38), it is evident that 



0 


(42) 
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This fact together with the transformations 

K = x 

ti = Py 

i = Pz 


p - 


m £ 


for < 1 reduces equation (37) to a Poisson equation like equation 
(l7b) , except that the independent variables are £ , q, and ( and the 
coefficient l/p^ appears on the right side of the equation. For heat 

addition of the type specified by equations (15) and (22) , cp(°) satis- 
fies Laplace's equation everywhere except precisely at x = 0. Thus, 

<p(o) 

must also satisfy the compatibility condition 


(4 0)+ - *« 0, "L> - 


The zero-order solution then is 
x,y,z system can be written as 


to equations (23) and in the 




On 


jmctx 


2icna3 


cos imtxz 


(43a) 


n==l 


,( 0 )- . - X 


On 


jcnox 


rt=l 


2rmap 


cos jmaz 


(43b) 


This solution is independent of y as expected, since to the zero or- 
der the effects of the body force are not included. Recall that the 
restriction m£ < 1 has already been imposed. Therefore, the factors 

3 and in the deno mina tors of equations (43) and others in this 

section should cause no concern. 

In order to solve the first-order system, substitute equations 
(37) into equation (40), integrate, and apply the boundary conditions 
to yield 


+ x l(x)g* 


(44) 
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to 

i- 


N 

O 


Since q* is independent of y, is also, and equation (39) "becomes 

V 2 ^ 1 ) - M 2 ^ 1 ) = 0 (45) 

tn XX 

Prom equations (41b) and (44) , the solution for <p(-0 involves the re- 
moval of the singularity at x = 0, just as in the small Mach number 
case previously discussed. Therefore, in a manner similar to the pre- 


vious method (eqs. (41b) and (44)), let 

<P^ = ? - M^yq^ 0 ) - x l(x)h(y,z) (46) 

where h is given by equation (24). Combining equations (41), (45), 
and (46) yields 

P \x + *yy + *zz = 5(x)j^yq* + P 2 h(y,z)J (47) 

p(—^y^z) = 0 (48a) 

<Py( x .>0,z) = f y (x,l,z) = 0 (48b) 

? z (x 5 y^0) » ? z (x,y,|) = 0 (48c) 


Applying equation (24) and the previously given transformation to the 
coordinate system to equation (47) yields 


OO CO 



C^cos 


JOUT] 

T 


COS 


TWJxt, 

6 


(49) 


where v* 2 denotes the Laplacian operator for the independent variables 
£, n, and £ and where the Pourier series representation 


oo 


z 

m=0 


G mn cos 


jtmT) 

P 



P 2 sinh 



rtna cosh 


mna 

2 
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with 


mn 


GM 2 

±1 oo 
2 


(m = 0) (50a) 


= 0 


4 «n 


M 


jt 2 jm 2 (m 2 + a 2 n 2 )J 


(m even) 
(m odd) 


(50b) 

(50c) 


is used in order that the boundary conditions on <p can be satisfied. 
The solution <j> once again satisfies Laplace's equation except pre- 
cisely at ? = 0 so that the following compatibility condition must 
also be satisfied: 


OO oo 


(*« - *4 =0 = 7 XX °® cos T ° os t 4 


rt=0 m=>0 


Hence, in the x,y,z system 


V* = IWCqo - 


T ~'\jm Z +(cm) 2 


P‘ 


m+n=l 


‘yjm 2 + (ctn) 2 


cos jtmy cos jtnozj 


(51) 


Thus, from equations (36), (43), (46), (24), and (5l) the function 
q> for large Froude numbers is to the first order in \ 




On T 


itnox 


n=l 


2jmop 


e P 


cos mnocz 


, x * 1 t x ) c 00 ( L y-V* Cma T S -Un2 + ( m )2\ 
, , •ST' On slnl1 ,mj '(y - I) 

(cos smy cos «nos) - x l(x) ^ ^ cos xnc 


itna jma 

n=0 cosh ~ 


( 52 ) 
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For zero Mach number, this solution reduces to that given for small Mach 
numbers by equation (28) . 


Result for uniform heat addition . - For uniform heat addition, as 
shown in figure 2(a), equation (22) becomes q*(z) = Qq with all other 

being zero. Then using equation (50) and applying L'Hospital's rule 
to equation (52) yield 


<p^= Qq! 


jM 


p* 


(l - Xy + |) + 3^ — -r — cos(2j + l)aty| 

V It 3 p (2j + l) 3 


Also, from equations (36), (42), and (44) 

XQqX l(x) 


* = 


P" 


(53) 


(54) 


From equations (12), (53), and (54) the velocity components are 


^ = 2° 
P 2 


CO 

1W (i + 


•x it(2J+l)x 

P 


U (2J+1)/ 


cos(2j + l)jty 


v* = 


2XQq 

* 2 P 


a> 

z 

>0 


ff(2J+l)x 

s_ P_ 

(23 + l) 5 


sin(2J + l)iry 


(55a) 


(55b) 


Except for the factor l/p in the coefficient and exponent of 
equation (55b), this equation and equation (30b) are identical. With 
increasing Mach number, then, and large Froude number, the total down- 
ward displacement of a streamline remains the same as in the sm a l l Mach 
number case; however, the axial distance over which the displacement 
occurs is decreased. 


The jump across x = 0 in the dimensionless increment of through 
velocity u is, from equation (55a), equal to Q q/P 2 J thus, the Mach 

number increases the jump. The series term in equation (55a), which, 
as has been pointed out, corresponds to a potential, also is increased 
with Mach number but dies out because of Mach number effects at shorter 
distances on either side of the plane x =» 0. 
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The component of u 


u =3 



represents a shear flow downstream of the heat addition and, again, is 
increased by Mach number. 


Result for cosine-wave heat addition . - In order to find the ef- 
fects of nonuniform heat addition in a plane transverse to the through 
flow, a cosine-wave heat-source distribution is assumed (fig. 2(b)), 
and q*(z) is given by equation (32). For this case, equation (52) 
becomes 


<P' 


,± „ § [ 

P 2 


1- X 


P 2 sihh 2jta 


fi - m£y) 


y " 2 


2jtax 


( y X) 


2jta cosh jta 


+ e P cos 2jtas /- 


CO 

X 


cos 2jtaz 


jtmx 


C _ ar 

mO + R 

^ — e p cos 
2jt ^ m 

m=l 


^ cos 2 itaz y . 

m=0 


"m2 


T~-\jm 2 +( 2a) 2 


~\Jm 2 + (2 a) 2 


cos nmy (56) 


and, from equations (36), (42), and (44), 


XQ x l(x) 

2 p2 


(l - p 2 cos 2itaz) 


-p 2jtOX 

M 2 e P 
00 

+ 4jtaP 


cos 2jtaz 


(57) 


The velocity components, determined from equations (12), (56), and (57), 
once again are composed of potential flows ■which vanish at distances of 
the order of H and L upstream and downstream of the heat-addition 
zone and the following components which persist for large x: 


u = 



v = - 



sinh 2jta 



2ita cosh na 


cos 2itaz 



cosh 2jta 



cosh Jta 


cos 2jtaz 


(58a) 


(58b) 
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w = - 


QXx 


sinh 2jta 


( y - t) ... 


cosh it a 


sin 2jtaz 


(58c) 


The v and w components are identical to those found in the previous 
section (eqs. (34h) and (34c)) and are independent of Mach number. The 
only effects of Mach number are apparent in equation (58a) , which de- 
fines a shear flow in this case which corresponds to that determined for 
uniform heat addition. The Mach number, however, increases the average 
(over the channel width) shear flow as compared with that for small Mach 
number (see eq. (34a)). The jump in axial velocity is also increased 
with Mach number. Equations (58) satisfy the simple continuity equation 
u^. + v y + w z = 0 and represent a vortex motion having components of 

vorticity 


& x = -jtQxaA sin 2jtaz 


(59a) 



2 , 


0 


— — 

jr - COS 2rtaz 


(59b) 

(59c) 


The x-component is identical to that given by equation (35a) for small 
Mach numbers , but the z -component is decreased by the Mach number. 
Therefore, although the cross-sectional streamline pattern is unchanged 
by Mach number (see fig. 5), the slope of the vortex filaments 9.^/0^ 


is infinite at the walls, as shown on the following sketch, instead of 



extending indefinitely downstream, as in figure 5. 
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Helicopter Rotor-Tip Ram-Jet Combustors 

The solutions of the previous sections are all linearly dependent 
on X. Therefore, the order of magnitude of X will be determined in 
a practical problem. For this determination, consideration is given 
to the rotor-tip combustors of helicopter ram jets. It must, of course, 
be realized that the configuration analyzed herein is only an idealiza- 
tion of an actual unit of this kind. Thus, 



where Urp is the tip speed, r is the rotor radius, and U_, of course, 

is the speed of the gas at which combustion is to take place. It is 
clear that the ratio H/r will be small in a practical case (H/r = l/l6 
would be a possible value), whereas Urp/U^ will be large, assuming that 

the flow at the inlet must be diffused to low velocity in the combustor 
(a diffusion ratio of Drp/CJ^ = 4 would be possible). Thus, values of 

X of unit order magnitude are not unreasonable. Actually, even if X 
is small, the analysis indicates that vortex motions can occur with 
nonuniform heat addition. 

For the case of uniform heat addition, the order of magnitude of 
the nonuniformity of the resulting flow can be estimated from the ratio 
of the velocity difference between the bottom and top of the channel D 
to the jump in velocity across the plane of heat addition AU- This ra- 
tio d/AH is just equal to X. For nonuniform heat addition, the ratio 
r/(2H + L)AU will yield an order of magnitude for the flow nonuniform- 
ity where r is the circulation in half the channel cross section and 
the bar over AH indicates that this is the average over z. Thus 

F x 

(2H + L)Z5J ~ X (1 + L/2H) 

Therefore, in both cases the nonuniformity in the flow evidently will 
be of the same order as X and hence can be appreciable as is demon- 
strated herein for helicopter ram- jet tip combustors. 


SUMMARY OF RESULTS 

An analysis of the steady, compressible, inviscid c h a nn el flow to 
which heat is a dde d at a cross-sectional plane and which is subject to 
a transverse body force showed how the flow depends on the Mach and 
Froude numbers and a dimensionless heat parameter. Solutions were 
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obtained for uniform or nonuni. form heat addition at a cross-sectional 
plane under the simplifying assumption of either smal 1 Mach number or 
of large Froude number. For small Mach numbers and uniform heat addi- 
tion, the flow was displaced in the direction of. the body force with 
the velocity profile being one of uniform shear. Nonuniform heat ad- 
dition in the cross-sectional plane in this case lead to the superpo- 
sition of strong vortex motions on the shear flow. 

The solutions for large Froude numbers yielded flows which were 
qualitatively the same as for small Mach numbers. For large Froude 
numbers, the shear component increased with Mach number, and the 
streamline displacement was more abrupt at higher Mach numbers. 

The flow .in a rotor-tip combustor of a helicopter ram jet wa s 
considered as an example, and the flow was shown to be appreciably 
affected by the action of a transverse body force. 


Lewis Flight Propulsion Laboratory 

National Advisory Committee for Aeronautics 
Cleveland, Ohio, November 10, 1955 
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APPENDIX - SYMBOLS 

The following symbols are used in "this report: 
constants defined Ly eqs. (25) 
constants defined "by eqs. (50) 
specific heat at constant pressure 
specific heat at constant volume 

velocity difference between bottom and top of channel for 
uniform flow 


Froude number , U t ,/-^/(-f 2 )H 
,th 


f i component of body force per unit mass, i = 1,2,3 
H channel height 

h function defined by eq. (18b) 

L channel width 

M Mach number, u/-y/r®T 

P pressure 

Q rate of heat addition by sources per uni t. volume 

On constants defined by eq. (22) 

Q constant defined by eq. (32) 

q,q* dimensionless heat parameters 

R gas constant 

r rotor radius 

T temperature 

U,U^,V,W velocity components, i = 1,2,3 

u,v,v dimensionless perturbation velocity components 
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X,X_£,Y,Z dimensionless coordinates, i =« 1,2,3 

x,y,z dimensionless coordinates 


O) 

CJl 

to 

w 


l(x) 


a 

P 

r 

r 

8(x) 


X 

p 

a 

T 

<P 

<P 

^i 

<p(i) 

X 

* 


unit step function, / 5(x)dx 

<J-ac 

x < 0 


length' ratio, H/L 


1 for x > 0, =» 0 


V 1 - 

circulation 

ratio of specific heats 


Dirac delta function, 1 ^ 


0 x ^ 0 


r 


x = 0 



6(x)dx 


1 


small parameter (see eq. ( 4 )) 
dimensionless parameter, l/lY^ a (-f 2 )H/u^ 
transformed dimensionless coordinates 
density 

dimensionless perturbation pressure 
dimensionless perturbation temperature 
potential function 

potential function defined by eq. (46 ) 

potential functions defined by eqs. (17), (18a), (19) 
(20) i = 1,2, 3, 4 


potential functions defined by eq. (36) 

dimensionless perturbation density 

density perturbation function defined by eq. (10) 


for 


and 
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*(i) density perturbation functions defined by eq. (36) 

£> x ,2y,£> z dimensionless vorticity components 
Subscripts : 

r reference condition 

T tip 

o° upstream condition 
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(a) Uniform heat addition; q* = Qq. (h) Nonuniform heat addition; 

q a 2 (1 - cos 2itaz) . 

u 

Figure 2. - CrooB-sectional pattern of heat addition at x =■ 0. Dark areas represent heat 
sources. 
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Figure 4. - Profile of increment in through velocity at large x for nonuniform heat 
addition. Mach number very much less than 1; length ratio, 1. 
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